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Abstract: 

We give a complete re-analysis of the leading twist quark-antiquark light-cone distribu- 
tion amplitudes of longitudinal and transverse p mesons. We derive Wandzura-Wilczek 
type relations between different distributions and update the coefficients in their confor- 
mal expansion using QCD sum rules including next-to-leading order radiative corrections. 
We find that the distribution amplitudes of quarks inside longitudinally and transversely 
polarized p mesons have a similar shape, which is in contradiction to previous analyses. 
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1 Introduction 



The theoretical interest in leading twist light-cone distribution amplitudes of hadrons is due 
to their role in the QCD description of hard exclusive processes [Jj . In terms of the Bethe- 
Salpeter wave functions these distributions are defined by keeping track of the momentum 
fraction x and integrating out the dependence on the transversal momentum k±: 

(f>(x) ~ / d 2 k ± 0(x,k ± ). (1.1) 

Jk\<p? 

They describe probability amplitudes to find the hadron in a state with minimum number 
of Fock constituents and at small transverse separation (which provides an ultraviolet (UV) 
cut-off). The dependence on the UV cut-off (scale) p is given by Brodsky-Lepage evolution 
equations and can be calculated in perturbative QCD, while the distribution amplitudes 
at a certain low scale provide the necessary non-perturbative input for a rigorous QCD 
treatment of exclusive reactions with large momentum transfer [0] . 

Their investigation has been the subject of numerous studies. Chernyak and Zhitnitsky 
(CZ) have developed an approach to study moments of light-cone distributions using QCD 
sum rules 0. Their main conclusion was j|, |5j that the pion and nucleon distribution 
amplitudes deviate strongly from the asymptotic distributions at large scales, which is a 
result still under debate. Another result ||, H was that the distribution amplitudes of 
longitudinally and transversely polarized p mesons deviate from their asymptotic distribu- 
tions in opposite directions: the longitudinal distribution is more wide while the transverse 
one is more narrow. In the further discussion, the pion and nucleon distributions received 
most attention. 

The present paper is devoted to the re-evaluation of p meson distributions along the 
lines of the approach of CZ and is mainly fuelled by newly emerged applications of light- 
cone distributions for the description of diffractive leptoproduction of vector mesons at 
HERA W\ and light-cone QCD sum rules for exclusive semileptonic B — > pev and radiative 
B — > p7 weak decays ||. The necessity of such an update is due to the following: 

First, the old calculations in |3], |f| have used a very low normalization scale p? ~ 
0.5 GeV 2 and a small value of the QCD coupling. Radiative corrections to the sum rules 
have been neglected. With the larger values of a s accepted nowadays the inclusion of the 
0(a s ) corections to the sum rules is mandatory. The corresponding calculation is a new 
theoretical result of this paper. 

Second, there is a controversy in the sign of the contribution of four-fermion operators 
to the sum rule for the transverse vector meson as given by CZ 0, [| , and later calculations 
||. This sign difference apparently remained unnoticed and has dramatic consequences for 
the shape of the distribution. 

Third, earlier studies have not spent due attention to distributions of transversely 
polarized quarks in longitudinally polarized mesons. As first noted in ||, to leading twist 
accuracy these distributions are given in terms of longitudinal quark spin distributions. We 
present a detailed derivation of the corresponding relations, the status of which is identical 



to that of the Wandzura-Wilczek relations [10] between the polarized nuleon structure 
functions gi(x,Q 2 ) and g2(x,Q 2 ). 

Our presentation is organized as follows. In Sec. 2 we collect relevant definitions and 
give basic formulas for the expansion of the distribution amplitudes in contributions of 
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conformal operators, which diagonalize the mixing matrix (Brodsky-Lepage kernels) to 
leading logarithmic accuracy. Section 3 is devoted to the analysis of QCD sum rules for 
the distributions in the transversely polarised p meson, while Sec. 4 contains the sum rules 
for the longitudinally polarized p meson. The final Sec. 5 contains a summary and some 
concluding remarks. We also include two appendices with the discussion of more technical 
issues. 



2 The p Meson Distribution Amplitudes 
2.1 Definitions 

We define the light-cone distributions as matrix elements of quark-antiquark non-local 
gauge invariant operators at light-like separations ||. For definiteness we consider the p + 
meson distributions; the difference to p° and u is just a trivial isospin factor in the overall 
normalization. The complete set of distributions to leading twist accuracy involves four 
wave functions 13b 



(0\u(0)v,»d(x)\p + (p,\)) 
(0\u(0) lp d(x)\p + (p,\)) 



(e ( %) 
(px) 



f p m p / due iup:c 0||(w,/i) 



(2.1) 



+ 



(A) 



X 



(px) 



f p m p l due-^gi\u,n), (2.2) 



(0\u(0) lpl5 d(x)\p + (p,\)} 
where the gauge factors 



W v p p x°f p m p f du e-^gf (it, p) , (2.3) 



Pexp 



ig \ da x^ A J ax) 
Jo 



are understood in between the quark fields. 

In the above definitions p p and are the momentum and the polarization vector of the 
p meson, respectively. The integration variable u corresponds to the momentum fraction 
carried by the quark. The normalization constants f p and f p (to be detailed later) are 
chosen in such a way that 

du f(u) = 1 

for all the four distributions f = 4>±, 4>\\, 9x\ 9± ■ The functions 4>±(u, p) and <p\\(u, p) give 
the leading twist distributions in the fraction of total momentum carried by the quark 
in transversely and longitudinally polarized mesons, respectively. The functions g±\u,p) 
and g±\u,p) describe transverse polarizations of quarks in the longitudinally polarized 
mesons and are to a large extent analogous to the spin structure function g 2 (x,Q 2 ) in 
polarized lepton-nucleon scattering. Similarly to the latter, they receive contributions of 
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both leading twist 2 and non- leading twist 3, and the twist 2 contributions are related to 
the longitudinal distribution 0y (u, y) by Wandzura-Wilczek |10| type relations: 



(v), twist 2 , 

9± I 



(a), twist 2 



(u,/i) = 2 



ul \ v h^ + u f 1 dv h^ 



(2.4) 



Here and below v = 1 — v etc. Eq. ( |2.4j ) is derived in App. A and presents one of our main 
results. 

The remaining twist 3 contributions to g±, g± can be written in terms of three-particle 
quark- ant iquark-gluon wave functions of transversely polarized vector mesons, cf. |3|, [Tl| , 
and will not be considered in this paper. From now on we will drop the superscript "twist 
2" , which is always implied. 

For some applications it is more convenient to rewrite (|2.2|) as 

(0|m(0) 7 ^(x)| P + (j9,A)) = p p (e^x)f p m p f due'^^u, y) 

+ e^f p m p fdue-tVgPfan), (2.5) 

J 



introducing a new distribution function 

1 



$nK/i) 



(2.6) 



Eq. (|2.6|) follows directly from ( |2.4|) and ( p75|) by integration by parts 



2.2 Conformal Expansion and Renormalization 



The separation between the quark and the antiquark in Eqs. ( p.l|) -( pT3|) is assumed to 
be light-like, i.e. x 2 = 0. Extracting the leading behaviour of the matrix elements on 
the light-cone one encounters UV divergences, whose regularization yields a non-trivial 
scale-dependence which can be described by renormalization group methods |@, |TJ. The 
conformal invariance of QCD at tree level implies that operators with different conformal 
spin do not mix with each other to leading logarithmic accuracy. For the leading twist 
distributions <p±(u,fi) and (f)\\(u, y) it follows that the coefficients a n of their expansion in 
Gegenbauer polynomials C%/ 2 (x) ]12| (that is in contributions of operators with definite 
conformal spin) are renormalized multiplicatively to that accuracy: 



(u, /i) = 6w(l — u) 



1+ £ aMC 3 J 2 (2u - 1) 



n=2,4,. 



a n (n) = a n (/i ) 



a s (/io) 



(7(n)-7(0))/(2A>) 



(2.7) 
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with (3o = 11 — (2/3) n/. The one- loop anomalous dimensions are [TJ| 

2 



n+l 



"5 I 1 - („ + !)"(„ + 2) + 4 g^' 



n+l 

7(n) = ^(l+4E 1 /j 

i=2 



(2.* 



The conformal expansion of the distributions g±\ gf 1 is more complicated and was 



derived in []8| using the approach of Refs. |TJ], IS]. We do not repeat the result in this 
paper, since to leading twist accuracy these distributions are not independent functions, 
but can be expressed in terms of <f>»(u,fi). One finds: 



g±\u,fi) = 6u(l - u) 



1+q4(/^2 /2 (0 + 



g { l \u, /*) = 2 (i + e 2 ) + A (ise 4 - ee - 1) + • • • , 



= -u{l-u)i 



1 + ^4(^(15^-11) + 



(2.9) 



Here and below we use the notation £ = 2u — 1 as shorthand. The leading contributions in 
Q2.9| ) agree with the "asymptotic distributions" that were derived in Ref. |TJ] by a different 
method, but erroneously identified as being of twist 3.[] 



2.3 Non-Perturbative Input 

The decay constants f p , and the coefficients a n in the Gegenbauer expansion ( |2.7|) are 
intrinsic hadronic quantities and must be determined either experimentally or by non- 
perturbative methods. In particular, the decay constant f p is measured [16|, [l7| : 

/p± = (195 ± 7) MeV, / p o = (216 ± 5) MeV. (2.10) 

For other quantities, most of the existing information comes from QCD sum rules. In 
what follows we summarize and update these calculations, taking into account radiative 
corrections and resolving some discrepancies in earlier studies. 



3 Transversely Polarized p Mesons 
3.1 The Tensor Coupling 

The normalization of the leading twist quark-antiquark distribution in the transversely 
polarized p meson is determined by the tensor coupling _/+, defined by 

(0\ua^d\p + (P, A)) = Ke^Pv - e^p^tf , (3.1) 

1 It is worthwhile to note that these leading terms correspond to the sum of contributions of leading 
and next-to-leading conformal spin, see || . 
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which can be estimated by studying the correlation function of two tensor currents within 
the framework of QCD sum rules ||18|| . We refer the reader to the reviews [[HJ and || 
for detailed explanations of the method; the latter reference deals specifically with the 
determination of distribution functions. A somewhat troublesome point in studying f~- is 
that the tensor current also couples to the positive parity J PC = l + ~ state b\ (1235) |jTT| : 



(0\^ lv d\bt(p i X))=ieJ^e^ Pfi fi 



(3.2) 

The correlation function of two tensor currents thus contains two Lorentz structures: 



n 



i Jd 4 ye iqv (0\T[u(y)a^d(y)d(0)a v ^u(0))\0) 



-z[(qx)(qpX v + q v Xy) - (qxfg^]I\. (q 2 



+ -^[(qx)(q^x u + q v x^) - {qxfg^ - q 2 x fl x l/ ]Tl + (q 2 



(3.3) 



To compactify the Lorentz structure we have contracted the correlation function in two 
indices by the light-like vector x^ ||. The Il ± (g 2 ) were calculated in || and correspond to 
intermediate states with positive (negative) parity, respectively: 



IF(g 2 



8vr 2 



q 2 In 



l + ^(ln 

37T v 



-<r 



7 

+ 3 



24q 



ra s qq) 2 



71 



9q 4 



4 

- ± 1 

9 



(3.4) 

where we used vacuum saturation for the contributions of four-fermion operators. 

The correlation function Il~(g 2 ) can be used to extract the value of f^~, see e.g. |§. 
Note, however, that it has a higher dimension than the correlation function of vector 



currents |18[], since in the latter case current conservation allows one to include one power 
of q 2 in the Lorentz structure. The higher dimension significantly reduces the accuracy 
of the sum rule, as it increases its sensitivity to higher resonances and the continuum. In 
addition, the sign of the four-quark contribution is reversed, which does not allow to get 
a stable sum rule for the p meson mass in this case, see 0. To overcome this difficulty, 
Chernyak and Zhitnitsky suggested to sum contributions of opposite parities. Since one 
has to assume 

IT + (0) + IT(0) = 

to avoid an unphysical singularity at q 2 = in Eq. ( |3.3|) , it is legitimate to write a dispersion 
relation for the structure 



n-( g 2 ) + n+(g 2 

q 2 



4vr 2 



In 



/i 2 



1 + 



3tt 



(in 



-q 2 7 
/i 2 + 3 



,a x ^ 2 1287T — _ 2 



12g 4 



tt G)+ 81q* 



. (3 - 5) 

Chernyak and Zhitnitsky speculated that the approximation of local duality for the 
continuum contributions may be satisfied with better accuracy in sum rules with summa- 
tion over different parity contributions, and noted that an additional advantage of using 



2 _B(1235) in old classification. 
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( p.5| ) is that contributions of particular four-fermion operators that are suspected to vio- 
late vacuum saturation cancel identically in this case. The price to pay is that the sum 
rule contains an additional contribution of the b\ (1235) meson; since its mass, however, is 
very close to the continuum threshold in the p meson channel, one may expect that this 
contamination has a minor effect. 

One can thus write down several different sum rules for f^, each of which has its own 
advantages and disadvantages, and their agreement indicates consistency of the approach. 
Using (|3.5|) one obtains 



(dse- s ' M 



4tt z 



1 + 



O'.s 



7T 



7 2 , s 
+ - In 



12M 2 [ n ' 81M 4 



9 3 p 

On the other hand, starting from the correlation functions IT T (g 2 ), one gets 



{yfcqq). (3-6) 



x 



1 + — 

7T 



7 2 , s 
- + - In — 

9 3 p 2 



1 ,a s „ 2 . 2087T ^_ 



™\2 



(3.7) 



7T 



7 2 , s 
+ - In — 
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80rr 



(V^qq) 2 , (3.8) 



where sg ~ 1.5 GeV 2 @ and s£ 



81M 2 

are the continuum thresholds in the p and 



2.3 GeV 2 

channels, respectively. The continuum threshold Sq for the "mixed parity" sum rule (|3.6j ) 
is discussed below. M 2 is the Borel-parameter. Note that the sign of the contribution of 
four-fermion operators in (|3.6| ) is opposite to the result given in || || . We have recalculated 
this contribution and confirm the sign as obtained in pj. 

In the numerical analysis we use a s (p = lGeV) = 0.56, i.e. = 0.4 GeV, corre- 
sponding to the world average a s {m z ) = 0.119 [JD]]. For the condensates we take the 
standard values WE 



.a 



-1 G ) = (0.012 ± 0.006) GeV 



7T 



^T s qq) 2 = 0.56 (-0.25 GeV) 6 . 



(3.9) 



The sum rules and the couplings are evaluated at p — 1 GeV. We have checked that 
changing the scale in the range p 2 = (1 — 2) GeV 2 does not have any noticeable effect, 
provided the extracted couplings are related by renormalization group scaling:^ 



/^lGeV) 



a s (lGeV) 
a s (p) 



4/27 
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Figure 1: (a) /^(l GeV) from Eq. ( |3.7|) as function of the Borel parameter M 2 for different 
values of the continuum threshold s . (b) The same for /^(l GeV) from Eq. Q3.8Q . 
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Figure 2: (a) f^(l GeV) from Eq. ( |3.6| ) as function of the Borel parameter M 2 for different 
values of the continuum threshold s . /^(lGeV) is put to 180 MeV. (b) The same with 



the bi contribution put to the continuum, i.e. fjj: = 0. 
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We start with the "pure parity" sum rules in Eqs. ( |3.7| ) and (|3.8|). The values of the 
couplings extracted from these sum rules are shown in Fig. |l|a and Fig. |l]b as functions 
of the Borel parameter for several choices of the continuum thresholds. Requiring best 
stability in the "working window" of the Borel parameter 1 < M 2 < 1.5 GeV 2 , we find 

ff = (160 ± 10) MeV, sg = 1.5GeV 2 , (3.10) 
f£ = 180 MeV, si 1 = 2.7 GeV 2 . (3.11) 



Note that s$ coincides with the value quoted in [18], while for b\ we get a somewhat larger 
value than Ref. 0. This difference, however, affects the coupling only very slightly: with 
Sq 1 = 2.3 GeV 2 we get = 170 MeV with a somewhat worse stability. Note also that it 
is difficult to specify more precisely the value of the continuum threshold s(J: the stability 
of the sum rule does not change much with Sq in the interval (1.3 -1.5) GeV 2 (although 
the value of /r~ does), which is precisely the disadvantage of having a sum rule of high 
dimension. 



Turning to the "mixed parity" sum rule ( |3.6|) we first note that the contribution of b\ 
is numerically supressed by the exponential factor exp[(m 2 l — m 2 )/M 2 ], so that a modest 
accuracy in is sufficient. Using the value in ( |3.11| ) as input and requiring best stability 



of the sum rule (|3.6| ) by varying M 2 and the continuum threshold (see Fig. |2|a), we get 

ff = (163 ± 5) MeV, s = 2.1GeV 2 . (3.12) 

The higher value of s (compared to Sq) is expected, since the part of the continuum 
contribution coming from b\ is taken into account explicitly on the left-hand side of the 
sum rule ( |3.6|) . 

On the other hand, since the b\ state is rather wide and its mass is very close to the 
continuum threshold in the pure 1 channel, it is natural to expect that an equally good 
fit to the sum rule can be obtained by ignoring this contribution on the left-hand side of 
( |3.6| ) and fitting the value of the continuum threshold to include it effectively. Remarkably, 
in this case we find a very similar value for the p coupling, see Fig. 

— (160 ± 15) MeV, Sq = (1.0 ± 0.2) GeV 2 . (3.13) 

Note that the sum rule now "wants" a much lower value of Sq. It is instructive to observe 
that the accuracy is now worse since the sum rule remains stable for a rather large interval 
of Sq- This is natural, since in this case we do not incorporate an additional information 
about the bi meson contribution. 

To summarize, we find that the positive parity b\ meson contributes significantly to 
the "mixed parity" sum rule, but it is not possible to separate this contribution from the 
continuum. In effect, the admixture of positive parity states can be described by lowering 
the duality interval for the p meson to 1 GeV. Our final result for the p meson tensor 
coupling is 

fp = (160 ± 10) MeV. (3.14) 



3 In this paper we stay consistently within 0(a s ) accuracy and do not attempt a renormalization group 
improvement of sum rules (see e.g. pQ]). 
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This value is by about 20% lower than CZ's result || and agrees surprisingly well with an 
old SU(6) symmetry relation, ff = (/ w + f p )/2 « 0.17 GeV [||. 

As discussed in [|, [J, an alternative method to determine /r - could be to consider the 
correlation function of the tensor with the vector current, which is not contaminated by 
positive parity states: 



dV iw {0\T[u(y)^d(y)d(0)a aP u(0)]\0) = [g^qp - g Ptl q a ] X (q 2 



(3.15) 



The correlation function was calculated in [22, [TT| and readsf] 



x(q 2 



2(qq) 



^(21n.^ 



m„ 1 
+ 3? + °'? 



(3.16) 



Here ml = (qgcrGq) / (qq). Note that the perturbative contribution vanishes to all orders 
and that the dimension seven operator qG 2 q has zero coefficient at tree level [22]. The 
corresponding sum rule reads (cf. P, |3|, p2|): 



e- m t /M2 fh») f P = ~2(qq) 




1 WIq 

3 



+ 



(9 2 S G 2 ) 
M 4 



(3.17) 



and yields « 200 MeV as illustrated in Fig. |. Here we use f p = 205 MeV, (qq)(l GeV) = 
(—0.25 GeV) 3 and ml = 0.65 GeV 2 at the scale 1 GeV. The accuracy of this sum rule is, 
however, not competitive to the ones above: the uncertainty in the quark condensate alone 
gives a 10% error; in addition, the study in [F2^] indicates possible large contributions of 
excited states to this sum rule, e.g. from p'(1600). Its significance is, however, that it allows 
to determine the relative sign of f p and f p , which proves to be positive. 



3.2 Deviations from the Asymptotic Form 

The deviation of the distribution function from its asymptotic form <p±(u) ~ u(l — u) is 
quantified by the coefficients a n in the expansion (|2.7[) . Since the corresponding anomalous 
dimensions are ordered with n, one can expect that, at least for large scales /x, only a few 
first terms are important. The QCD sum rule approach can be used to estimate a^- The 
traditional procedure developed by Chernyak and Zhitnitsky is to write down the sum rule 
for the second moment of the wave function, which is related to by simple algebra: 

r 1 1 12 

/ du{2u-l) 2 ( f )± (u, f i) = - + — 4Gu). (3.18) 
Jo 5 35 

The corresponding sum rule is obtained from the correlation function of the tensor current 
with the similar operator with two extra covariant derivatives u(y)a p ^x^(i D -x) 2 d(y). 
We find it more appropriate to consider the sum rules directly for the coefficients in the 

4 The radiative correction to the quark condensate contribution is a new result. 
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Figure 3: /^(lGeV) from Eq. ( [3.17 ) as function of the Borel parameter M 2 for sq 
1.5 GeV 2 . 



expansion in Gegenbauer polynomials, which in general correspond to correlation functions 
of the tensor current with the conformal operators 



n T ± (n \y) =i n (d.) n 



u(y)^.c 3 J 2 ( g— & ] d(y) 

D. + D. 



where the dots stand for the projection on the light-like vector x^ and d is the total 
derivative. Note that one of the indices of the a matrix is projected onto x M , while the 
other one has to be taken transverse to the (x, q) plane, where q is the p meson momentum, 
see 



15 for more details. 



As a general property of conformal operators [^3J the tree-level perturbative contribu- 
tion to the corresponding correlation function vanishes (for n ^ 0) and the perturbative 
contribution to the corresponding sum rule starts with order 0(a s ). As a result, these 
sum rules are necessarily less stable than the sum rules for moments, and their accuracy is 
seemingly worse. The better accuracy of the sum rules for moments is, however, completely 
illusory since in this case the major contribution comes from the trivial first term in ( |3.18| ), 
corresponding to the asymptotic distribution function, and the contribution of interest is 
numerically suppressed. Since we should not expect good stability for the sum rule for a2, 
we evaluate this sum rule using precisely the same values of the continuum threshold and 
the same "window" of the Borel parameter as in the sum rules for fr~. The instability of 
the sum rule then gives an estimate of the accuracy of the result .Q 

It is important to note that the necessity to separate the contribution of leading twist 
does not allow for the separation of contributions of opposite parity in the diagonal sum 
rules. Indeed, one may try to start from the correlation function like the one in (|3.3j) 



with two open Lorentz indices (and with the substitution of one of the tensor currents 

5 It has become a common practice to choose different values of the continuum threshold in the sum rules 
for different moments. To our point of view, the higher fitted values of so for higher moments n = 2, 4, . . . 
generally reflect the increase of the overall mass scale in the correlation function, due to the increasing 
contribution of higher resonances. This rise has nothing to do with the change of the interval of duality 
for the p meson contribution of interest, which is in fact more likely to decrease. 
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by ft^^(y)), and try to isolate the negative parity contribution by taking the projection 
q fJj q u H^j) n ^ = (qx) n+2 Tl~ T ( n \q 2 ). However, the same projection for the defining Eq. ( O 



vanishes identically since to leading twist accuracy one must put contributions of order 
q 2 = m 2 to zero. Thus, this projection is in fact saturated by higher twist contributions 
and is irrelevant for our analysis. Therefore, one cannot get rid of the contribution of states 



with positive parity and a more convenient correlation function to consider is: |LI 



i Jd 4 ye tqy (0\T[u(y)a^d(y)d(0)a^x ( (i D -x) n u(0)}\0) = -2(gx) n+2 n T(n) (g 2 ). (3.19) 

It is easy to check that the trace over Lorentz indices picks up the required transverse 
components. 

The complete results for the sum rules for the coefficients in the Gegenbauer expansion 
for arbitrary n are given in App. B. Note that in this case the mass scale in the correlation 
functions rises as M 2 ~ n 2 for large n as compared to the increase M 2 ~ n for the moments. 
This makes the sum rule approach essentially useless for the evaluation of a n with n > 2. 
For the particular case n = 2 we get, using the correlation function ( |3.19|) : 

-mj/M 



* IM \ft;?(p)^-^(p) +h meson = 



1 ° M «*[! " «-""] ■ U < <? > + £k < ^ > 2 . (3.20) 



2tt 2 7T L 1 5 3M 2 x 7T ' 9M 4 

This sum rule is equivalent to the sum rule for the second moment considered in || 



e- m ' /Ma (f^) 2 (u) Cdu (2u - l) 2 (P ± (u,fx) + b x meson 
p Jo 



d^~ S/M 1 + -T 77 + 21n-^ ) \ + ^tj, G 2 ) + — -A^qq) 2 ■ (3.21) 



20vr 2 y 1 7T V 15 /i 2 /J 36M 2 N 7T ' 81M 4 



provided one takes the same value of the continuum threshold as in the sum rule for the 
tensor coupling ( ^.6[ ). Note that the sign of the contribution of the four-quark condensate 
is opposite to the result of 0, [3|].[] 

The value of that follows from the sum rule ( |3.2U| ) is plotted as a function of the 



Borel parameter in Fig. f|. Note that we do not have an independent estimate for the 
contribution of the b\ meson in this case, so we neglect it and take a low value for the 
continuum threshold, Sq = 1 GeV 2 , on the right-hand side. From this we get as our final 
result 

40 = 1 GeV) = 0.2 ±0.1. (3.22) 

This has to be compared with a^~(/i = 1 GeV) = —0.17 from |J [|; the difference in sign is 
mainly due to the opposite sign in the contribution of the four-fermion operators in || [| . 

We have investigated whether adding the b\ contribution as a free parameter and re- 
quiring best stability in the range 1 < M 2 < 1.5 GeV 2 could change the result. We have 



6 We have recalculated this contribution and get the opposite sign for all moments, see App. B. For the 
case n — our result agrees with M. 
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M 2 [GeV] 2 

Figure 4: a^-(l GeV) from Eq. (|3.20|) as function of the Borel parameter M 2 for sq — 1 GeV 2 



also tried to follow the standard procedure to use the sum rule ( 13.21] ) with s fitted to get 
best stability. In both fits the value of tends to increase by some (30-50)%, but we do 
not find this evidence significant enough to influence our estimate. 

To avoid an admixture of positive parity states, one can consider, instead of ( |3.19| ), the 
correlation function 



Wye™ (0\T[u(y) ltl d(y)d(0)a a pX P (i D ^)>(0)]|0> = [g a ,{qx) - x,q a ] {qx) n X ^\q 2 ). 

(3.23) 

The results for x {q 2 ) are available from [2|J]. The corresponding sum rule for reads 



e- m * /M2 fMf P a^) = -^(qq) 
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, 29 a s I M 2 7ds s/M 2 

1 + — — In — - 7E + const. - / — e~ s/M 



- 2 



+ 



18 7T 

85 (g 2 G 2 ) 



M 2 216 M 4 



(3.24) 



where we used vacuum saturation for the contribution of dimension seven. The constant in 
the radiative correction to the quark condensate contribution is not calculated yet. Unfor- 
tunately, due to the large coefficient in front of the contribution of the mixed condensate, 
its contribution almost identically cancels the leading quark condensate contribution, and 
the answer depends crucially on the contribution of dimension seven, which is poorly known 
(it is suspected that vacuum saturation is strongly violated in this case). Thus, from this 
sum rule one can only get a rough estimate \a^ \ < 0.5. 



4 Longitudinally Polarized p Mesons 

Since the decay constant f p is measured experimentally (we use the average value f p = 
(205 ± 10) MeV in the numerical analysis), we only need an estimate of the coefficient a| 
describing the deviation of the distribution 0m from its asymptotic form. The corresponding 
QCD sum rule calculation has been done by Chernyak and Zhitnitsky in Ref. ||. We update 
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0.6 0.8 1 1.2 1.4 1.6 1.8 2 



M 2 [GeV] 2 

Figure 5: c4(lGeV) from Eq. ( |4.1|) as function of the Borel parameter M 2 for s = 
1.5 GeV 2 . 



this calculation by taking into account radiative 0(a s ) corrections and using an up-to-date 
value of the strong coupling that is slightly larger than the value used in Ref. @|. The 
radiative corrections can be extracted from a paper by Gorskii [25|, where he calculated 
the correlation function of two axial vector currents (with extra derivatives), which in 
perturbation theory and for massless quarks coincides with the vector correlation function. 
The complete results for arbitrary moments are given in App. B. For n = 2 we get the sum 
rule 

e-«"*fl % 4(ri = ± ^ M*[l-e-o/«° ]+ ^ ( ±<?) + ** < vW , (4.D 
which is equivalent to the sum rule for the second moment considered in [|J: 



e~ m V M2 f p jUii(2tt-l)Vo(«,/0 = 

20tt 2 I + 3 n ) y ) + 12M 2 { n U > + 81M 4 Wasqq) " 1 j 

With the same input parameters as in Sec. 3, the numerical analysis yields (see Fig. ^): 

af(/i= 1 GeV) = 0.18 ±0.10. (4.3) 

This is in perfect agreement with the original estimate a|(/i = 1.1 GeV) ~ 0.18 H. It 
also coincides within the errors with our result for a^~, Eq. (|3.22|) , which means that the 
distribution amplitudes <ft\\ and <f)± are similar. 



5 Summary and Conclusions 

Extending earlier studies 0, || 0], we have performed a re-analysis of p meson quark- 
antiquark light-cone distribution amplitudes of leading twist. In general, their complete 
set consists of four independent functions, but we have shown that to our (twist 2) accuracy 
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Figure 6: Final results for the wave functions <j)\\ (a) and <p± (b) at /x = 1 GeV (solid lines). 
Long dashes: asymptotic wave functions, short dashes: 0j_ according to CZ B. 



the distributions of transversely polarized quarks in the longitudinally polarized p mesons 
can be related to the distributions of longitudinally polarized quarks. The theoretical 
status of these relations is identical to the status of the Wandzura-Wilczek relation [10] 
between the polarized structure functions of the nucleon, gi(x, Q 2 ) and g2(x,Q 2 ). 

We have given a detailed re-analysis of the QCD sum rules for the first two moments of 
the distribution amplitudes, complementing existing sum rules by the calculation of radia- 
tive corrections. Our final results for the distribution amplitudes of quarks in longitudinally 
polarized and transversely polarized p mesons are shown in Fig. |6|a and b, respectively. The 
solid curves correspond to the distributions calculated using the parameters specified in 
the text, the dashed curves show the asymptotic distributions. 

We deviate from the results of || |J mainly in the shape of the distribution amplitude for 
the transversely polarized p meson, see the short-dashed curve in Fig. |b. We find that the 
distributions in longitudinally and transversely polarized p mesons coincide to our accuracy, 
whereas in || U a significant difference has been claimed. This contradiction is largely due 
to an opposite sign of the contribution of four-fermion operators in the corresponding sum 
rule. We note that the sign as given in || Q also contradicts an independent calculation 
in Ref. ||, which apparently remained unnoticed. One more consequence of this sign 
difference is that our result for the tensor coupling is 20% lower than in || § . 

A discussion of phenomenological consequences of our results goes beyond the tasks of 
this paper. Since in hard exclusive processes one typically deals with integrals over quark 
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distributions of type 



Jo 



du 



u(l — u) 



6f p [l + a 2 (Q) 



the change in shape of the transverse p distribution suggested by the results of this paper 
may increase the rate of the production of transversely polarized p mesons by a factor two. 
The consequences for exclusive semileptonic and radiative B decays will be considered in 
a separate publication [26 . 
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Appendix A: Transverse Spin Distributions 

The derivation of relations between longitudinal and transverse quark spin distributions 
in the longitudinally polarized p meson is in principle straightforward and can be done 
similarly to the classical Wandzura-Wilczek analysis for polarized leptoproduction |10 . 
The major difference is that one must include operators with total derivatives and that 
higher twist operators corresponding to total derivatives of lower twist operators cannot 
be neglected. 

It is convenient to consider the relevant non-local operator at symmetric quark- ant iquark 
separations: 



u(-x)^d{x) = x" 1 ■ ■ ■ x^—u{0) D m ■ ■ ■ J A ) ■ 



nl 



(A.l) 



Since x 2 = the arising local operators are traceless (contractions of the type vanish 
by the equations of motion), but not fully symmetric in Lorentz indices because of the 
distinguished index p. Therefore, they contain a mixture of contributions of twist 2 and 
twist 3, which have to be separated: 



u(-x)jA x ) = u{-x)^^d(x) 



twist 2 



+ 



u{-x)^^d(x) 



twist 3 



(A.2) 



The leading twist 2 contribution by definition contains contributions of symmetrized op- 
erators: 



w(-x)7 A1 d(x) 



twist 2 



n=0 



nl 



-u(0)< 



n + 1 



D 



/ii 



n 



Be.-, 7 fc }d(0) ■ (A.3) 



Fortunately, the sum can be re-expressed in terms of a non-local operator [27 

u{-x)^^d(x) 



twist 2 







1 d 
dv — — u(— vx) $d(vx) 



dx 



(A.4) 



A' 
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which is easily verified by expanding. An identical expression is valid for the non-local 
operator with an additional 75 in between the quarks. 

Using the equations of motion, the difference u(— x^^d^x) — [u(— x)j fl d(x)]t w i B t2 can 
be written in terms of operators containing total derivatives and quark-antiquark-gluon 
operators, see P7| . Neglecting quark masses, one finds: 

w(-x)7 jtt d(x) 



twist 3 



1 ru r _ -| 

g s I du dvu(-ux) uG llu (vx)x v ^75 — ivG llu (vx)x v d(ux) 



u(-xh^ 5 d(x) 



twist 3 



+ ie ^ / udux v d a u^—ux^^diux) 



rl pu r ~ 1 

—g s du dvu(-ux) uG pv {v x)x v $ — iv G^iv x)x u ^75 d(ux) 

JO J—u L ' - 



+ ie 1,01/3 / udux v d a u(—ux)^pd(ux) 



(A.5) 



where G pu is the gluon field strength, G^ = (l/2)e llva pG a(i ', and d a is the derivative over 
the total translation: 



u(—ux)^pd{ux) 



d 



u(—ux + y)j^d(ux + y) 



(A.6) 



y^O 



Note that ( A.4 ) and ( A.5 ) are exact operator relations. Taking the matrix element 
between the vacuum and the p meson state, we get 



u(~x)-f ll d(x)] tvf . st2 \p + (p,X)} 



1 d 

dv — — (0|w(— vx) $d(vx)\p + (p, A)) 
ox a 



1 , d 



dv ( e W x )f p7n I due-^^u) 



dx 



ei X) f P m p 



Jo 



dv / due-^ vpx c 



>\\(u) - ip^e^x 



u 



V P 



f p m p [du e-^>|| (u) + f e W - pj£$) f P m P i! du [ dv e ~^ll («) > ( A - 7 ) 



(px) 



Jo 



where £ = 2w — 1 and to arrive at the last line we have used 

( l dvv f l du£e~^ x (j)\\{u) = — C dvv f du-^-e-^c 
Jo Jo px Jo Jo Ov 



U I 



px 
i r 1 

— / du<p\\(u) 
px Jo 



iipx - / dve-^ vpx 



Note that the matrix element of the twist 2 operator produces both Lorentz structures, 
and hence g±\u, p) is nonzero to this accuracy. 

Specific for the kinematics in exclusive processes is the generation of an additional twist 
2 contribution by twist 3 operators proportional to the total derivative d a — > —ip a , which 
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would vanish in deep inelastic scattering. Taking the matrix element for the twist 3 operator 
in the first of Eqs. ( |A.5|) and neglecting three-particle quark-antiquark-gluon distributions 
of twist 3 [11] we get 

(0| u(-x)-ff,d(x) 



twist 3 



P + (P,A)) 



) f p m p C v 2 dv f 1 du e-***g<f> (u, p) . (A.9) 



Since, on the other hand, 



(px) 



(e< A >x) 



\u(-x)^ p d(x)\p + (p, A)) = p p — — ^f p m p due iipx (j)\\(u, y) 



(px) 



/ p m p / due*>°gP(u t ri, (A.10) 



we obtain relations between gf^(u,p), g±(u,[/,) and 0||(w,/i) by comparing the Lorentz 
structures. At this stage it is convenient to introduce the moments 



Ml = j o duCH K = / duZ n g^\u,n), Ml 

Expanding ( |A.7| ), ( |A.9| ), ( |A.10| ) in powers of (px), we get 

AC 



dur^W)- (A.n) 



1 run — 1 ) , 1 ,ii 

AC_, + tM\ 



2 n + 1 "~ n ~ z ' n + 1" n 
Similar manipulations with the axial-vector operator (|2.3j ) produce one more relation 

1 



(A.12) 



1 

-M a = 

2 n n+2 n 



-Ml. 



(A.13) 

Note that the contribution of the leading twist operator (0| [u(— x) r y Pj j5d(x)} tw ist2\p + (p, A)) 
vanishes identically in this case, and the answer is generated entirely by twist 3 operators, 
which are reduced to total derivatives. 

Combining ( [A. 12] ) and ( A.13 ) we get a simple recurrence relation, 



{n + l)M: = {n-l)M:_ 2 + Ml 



(A.14) 



the solution of which yields the first relation in (2.4). The second one then follows from 
( |A.13| ) after some algebra. 



Appendix B: QCD Sum Rules for Arbitrary Moments 

In this appendix we collect some more definitions and give the sum rules for the Gegen- 
bauer moments a n of the longitudinal and transversal p meson distribution amplitudes for 
arbitrary n. 
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We first relate the a n to hadronic matrix elements of local operators. To leading loga- 
rithmic accuracy, the relevant multiplicatively renormalizable operators are: 

n 

n v(n) {y) = Cn>j {ixdr^u{y) jt (ix D) s d(y) , 

3=0 

n _ 

tf {n \y) = Y, c n,A^d) n - j u{y)a^x v {ix D) j d(y) , (B.l) 

3=0 

<—> > ■* 

where x p is a light-like vector, = (i/2)[y^y v and D p =d P - <9 M -2igA p [(y)X a /2. 

The c n ,k are the coefficients of the Gegenbauer polynomials such that C^J 2 (x) — ^Cn^x^. 
Expanding the defining relations for 0, Eq. fl2.1|) -( |E3"l) , around the light-cone, one finds 



(0\Q v ^{0)\p) = {px) n f P m p {ex) duC 3 J 2 (2u - 1) <f> {l {u, p) 



{px) n f p m p {ex) 3(n 2 | 2 ^| n 3 | 2) 4 (p) 



<0|fiJ<»>(0)|/>) = {pxTif^{e,{px)-p,{ex)) Jdu C 3 J 2 (2u - 1) </> ± (u, p) 



= (pxTif^ipx) -^W)^^«nM. (B.2) 

The QCD sum rules |18| are obtained by matching the representation in terms of hadronic 
states to the operator product expansion in the Euclidian region for the correlation func- 
tions 

(gx) n+2 n y(n V) = i U A ye iqv (0|TO y(n) (y)ft ty(0) (0)>, 



-2(gx) n+2 n T ( n V) = i Jd*ye iqy (0\m T ^ n \y)^ T ^(0)}. (B.3) 

Note that the contraction over p in the second relation automatically projects onto the 
transverse component Q± \ which is a conformal invariant operator, whereas f^ n ) is not. 
We find the following sum rules for (for even n): 

KJp ) 2(2n + 3) n w 



1 

M 2 [l - e~ So/M2 ] JduuuCl /2 (2u - 1) (\nu + lnu + In 2 ~ 




1 a M M^U _.-VM 2 
2n 2 ix 



+ ^(^ G2 )( n2 + 3n - 2 ) + ^AV^m)Hn-l)(n + l)(n + 2)(n + A).(BA) 
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The radiative correction in ( |B.4j ) is a new result. Similarly, we obtain for a\: 



2 /M 2 2 3(w+ l)(w + 2) J, 
Jp 2{2n + 3) nW 



+ ^4 ( y/wi ) 2 (» + + 2 )(^ + 3n - 7) . (B.5) 

In this case a compact answer for the radiative corrections as in ( p.4j) is not available, 
but the rj[ are related to the radiative corrections to the axial vector correlation function 
(with extra derivatives) and for arbitrary n can be expressed in terms of the coefficients 
A' k calculated in fl25| : 



(* + !)(* + 3) 



r l = E C ^ Ti , m,^ • ( B - 6 ) 



In particular 



A) = l, ^2 = 3, 4 = 27' A6 = l3^ : 



1 11 83 



' — o ' 1 2 — o ' ' 4 — C ) '6 



3' z 3' "6' u 810 

For completeness and for comparison with || , we also give the sum rules for the moments 
(S n )=fdu£ n (P(u,ti): 



O S ° 1 ( ( 2 

(/^WDiM e~™V M2 = JL J ds jdue-»' M2 uu(2u - 1)" 1 + ^ 6 - ?L + 21n ± 

71 I ^ \ A* 

,2 M M ra — 1 1 , a s „ 2 , 647r — 2 

+ ln« + In + In -)} + __(_ G» ) + _ („ - 1) < ^ )>, (B.7) 



Note the difference in sign in the last term in Eq. ([B.7]) with respect to (4.25) in Ref. 0. 
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